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Fig.  1.  A  traditional  active  radar  system  with  one  target. 

I.  Introduction 

Compressive  sensing  (CS)  has  received  considerable  attention.  CS  has  shown  great  potential  in 
diverse  fields,  for  example,  high-resolution  image  processing,  wireless  communication,  and  radar 
systems  [l]-[6].  The  central  goal  of  CS  is  to  reconstruct  a  signal  using  very  few  measurements 
[7].  This  report  starts  with  a  brief  review  of  target  detection,  communication  and  CS  theory  for 
different  radar  and  wireless  communication  signals.  And  then  the  report  focuses  on  different  CS 
applications  and  perfect  signal  reconstruction  conditions. 

II.  Adaptive  Compressed  Sensing  Radar  Settings 
A.  Preliminary  on  Traditional  Radar 

We  begin  with  an  introduction  to  a  traditional  radar  system.  A  simplified  radar  system  is  given 
in  Fig.  1.  Such  a  system  is  one-dimensional,  monostotic,  single-pulse,  and  active  radar.  Suppose 
a  target  located  at  range  x  is  traveling  with  constant  velocity  v  and  has  reflection  coefficient  sxv. 
After  the  transmission  of  signal  f(t),  the  receiver  observes  the  reflected  signal  as 

r(t)  =  sxvf(t  -  rx)  exp(2yi -juxvt)  (1) 

where  tx  is  the  round  trip  time  of  the  signal,  j  =  \f—l,  and  ooxv  is  the  Doppler  shift.  The  target 
information  (x,  v )  can  be  readily  obtained  from  (rx,u>xv),  the  pair  of  ( time  delay ,  Doppler  shift), 
of  /(•).  Here,  an  ideal  case  is  considered  without  noise. 
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Using  a  matched  filter,  the  reflected  signal  r  is  correlated  with  a  time-frequency  shifted  version 


of  the  transmitted  signal  through 


(2) 


Jr 


where  Arf(-,  •)  is  so  called  the  cross-ambiguity  function  of  two  functions  /  and  r,  and  the  super¬ 
script  *  denotes  the  conjugate  operation.  From  (2),  we  can  see  the  ambiguity  surface  of  /  centered 
at  the  target’s  location  and  scaled  by  its  reflection  coefficient  |  | .  It  is  straightforward  to  extend 

(2)  to  include  multiple  targets.  However,  when  two  or  more  targets  are  too  close  to  each  other 
and  have  overlapping  ambiguity  functions,  there  are  some  uncertainties  on  the  exact  locations  of 
targets.  Therefore,  the  range- velocity  resolution  between  targets  of  classical  radar  is  limited  by  the 
radar  uncertainty  principle. 

B.  Fundamentals  on  CS  and  CS-Based  Radar  ( CSR) 

Under  appropriate  conditions,  CSR  can  beat  the  traditional  radar.  We  now  consider  I\  targets 
with  unknown  range-velocities  and  corresponding  reflection  coefficients  in  an  N'  x  N  grid  plane. 
If  the  number  of  targets  K  <C  N'N  and  the  transmitted  signals  are  incoherent,  we  can  use  the  CS 
technique  to  recover  such  a  sparse  target  scene. 

A  CSR  has  the  following  features: 

1)  Eliminating  the  need  of  matched  filter  at  the  receiver; 

2)  Requiring  low  sampling  bandwidth  for  low  information  rate  (other  than  the  high  Nyquist  rate); 

3)  Providing  higher  resolution  compared  to  the  traditional  radar. 

Considering  a  target  scene  in  an  L  x  M  range-Doppler  plane  in  Fig.  2,  we  define  the  time  delay 
and  Doppler  shift  matrices,  respectively,  as 


^  0  IxN  ^ 


(3) 
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Fig.  2.  The  range-Doppler  plane  of  the  target  scene  that  CSR  measures  is  discretized  into  an  L  x  M  grid.  The  target  located  on 
the  Ith  range  and  the  rath  Doppler  bin  is  labeled. 
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0  *  *  *  0  1  / 

where  L  and  M  denote  the  numbers  of  range  and  Doppler  bins  that  CSR  measures,  respectively, 
N  is  the  length  of  transmission  waveform  vector,  and  %  =  exp(27rj/M)  is  the  root  of  unity. 

A  basis  dictionary  H  =  ( | h0,o I  |ho,i|  •  •  •  |hL_i;M-i|)  can  be  created  by  concatenating  the  received 
signals  where  h/  m  =  Fm  •  T  •  s.  Here,  s  denotes  the  transmission  signal  of  length  N.  H  contains 
all  the  possible  signal  reflected  from  the  target  in  any  grid  of  the  range-Doppler  plane  and  is  named 
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as  sparse  representation  matrix  in  the  CS  theory. 

Define  9^m  as  the  complex  scattering  coefficient  in  the  (L  m) th  grid  of  the  target  scene,  and  the 
corresponding  coefficient  vector  is  denoted  as 

(  #0,0  #0,1  •  •  '  #L— 1,M— 1  )  (5) 

where  the  superscript T  denotes  the  transpose  operation.  The  task  of  a  CSR  is  to  discover/identify 
the  coefficients  in  the  range-Doppler  plane.  The  received  signal  reflected  from  the  target  scene  can 
be  expressed  as 

L—l  M—l 

r  =  +  n  =  HO  +  n  (6) 

1=0  m= 0 

where  n  is  the  additive  Gaussian  white  noise  (AWGN)  vector  of  length  L.  The  compressed  signal 
y  can  be  written  as 

y  =  $r  =  +  (7) 

T 

where  $  is  the  P  x  L  measurement  matrix. 

In  adaptive  CSR,  both  transmission  waveforms  s  and  measurement  matrix  are  treated  as 
variables  which  will  be  optimized  and  updated  based  on  0,  the  estimation  of  6.  Note  that  T  =  d>H 
is  the  equivalent  dictionary  for  this  proposed  adaptive  CSR.  Complex  Gram  matrix  G  is  applied  to 
represent  the  mutual  coherence  and  it  is  defined  as  [8] 


g  =  rHr 


(8) 


where  the  superscript  H  denotes  the  conjugate  transpose  operation. 

Variables  s  and  are  optimized  to  make  the  resulting  Gram  matrix  approximate  a  diagonal 
matrix  (the  ideal  case),  i.e.,  G  =  Diag(g0,o>  9o,i,  •  •  •  ,  Algorithms  are  proposed  for 

the  optimization  by  fixing  s  and/or  <E>.  With  the  changing  of  T,  the  target  is  to  makes  the  realistic 
Gram  matrix  G  as  close  to  the  ideal  diagonal  G  as  possible.  This  can  increase  the  resolution  of 
CSR  significantly.  In  Fig.  3  on  the  next  page,  a  work  flow  chart  is  given  where  G  is  the  estimation 
of  G. 
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Fig.  3.  Work  flow  chart  for  this  adaptive  CSR  system. 


III.  Analog-to-Information  in  CS 


In  this  section,  we  first  generalize  CS  system  for  continues-time  (CT)  signals,  and  then  we 
present  some  more  possible  applications  in  the  field  of  CS. 

A.  Traditional  Form  of  CS 

We  describe  the  traditional  form  of  a  CS  system  in  discrete-time  (DT)  domain  in  this  subsection. 
CS  deals  with  the  problem  of  acquiring  an  N  x  1  discrete  time  signal  vector  f .  To  proceed,  f  should 
be  K- sparse  in  some  N  x  N  sparsity  basis  matrix  in  which  each  column  is  a  basis  vector  xpi, 
i.e.,  ^  =  [t/q,  •  •  •  ,  i/^,  •  •  •  ,  If  we  represent 

f  =  ^x  (9) 

where  f  is  A'-sparse,  meaning  that  only  K  <C  N  expansion  coefficients  in  x  are  non-zero.  With  <b 
denoting  an  M  x  N  sensing  matrix,  we  have 

y  =  $f  (10) 

which  represents  the  iV-dimensional  input  signal  encoded  into  an  M-dimensional  set  of  measure¬ 
ments  and  M  <C  N.  Using  (9)  and  (10),  y  can  be  rewritten  as 

y  =  #f  =  #\j>  x  (11) 

u=#^ 

where  y  is  an  M  x  1  vector,  $  is  an  M  x  N  matrix,  ^  is  a  N  x  N  matrix,  and  x  is  a  N  x  1  vector. 

In  the  present  form,  CS  is  only  applicable  to  DT  signals.  Below  we  extend  the  framework  to  CT 

signals.  In  the  literature,  there  exist  a  few  kinds  such  conversion  architectures.  We  will  present  an 
architecture  of  such  conversion  from  CT  signal  f(t )  to  DT  signal  f  in  the  next  section. 

B.  Analog-to-information:  CS  from  CT  to  DT 

We  now  generalized  CS  theory  from  the  DT  to  CT  domain.  To  develop  CT  domain  CS,  one  will 
need  to  define  analog  signal  models  for  sparse  signals  and  then  construct  an  analog  system  which 
is  CS -compatible. 


7 


Pc,m(t) 


Fig.  4.  An  analog-to-information  CS  architecture. 


Assuming  the  analog  signal  /(f)  can  be  represented  through  a  finite  number  of  parameters  in 
some  continues  basis  ^(f),  i  =  1,  •  •  •  ,  TV,  then  /(f)  can  be  decomposed  as 

N 

fit )  =  Xi  G  R.  (12) 

1=1 

In  the  case  that  there  are  a  small  number  of  nonzero  entries  in  x  =  [sq,  x2,  ■  •  •  .  a-'Jv]T,  again  we 
can  say  that  the  signal  /(f)  is  sparse. 

The  next  question:  how  one  can  convert  CT  signal  /(f)  to  DT  measurements  ym  (m  =  1,  •  •  •  ,  M)? 
Figure  4  gives  an  architecture  to  solve  above  issue.  In  Fig.  4,  Pc>m(f)  denotes  the  mth  distinct 
chipping  sequence  which  is  a  pseudo-random  sequence  of  ±ls.  The  alternating  frequency  of  the 
chipping  sequence  must  between  values  at  or  faster  than  the  Nyquist  frequency  of  the  input  signal. 
That  is,  fcs  is  a  much  lower  (compared  to  the  Nyquist  rate  fs  =  2 /max)  sampling  rate  (e.g., 

fcs  =  fs/N). 

From  Fig.  4,  ym  can  be  obtained  as  a  result  of  modulation,  convolution  and  sampling  (at  rate  of 
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fm )  operations.  Therefore,  we  can  write  ym  in  the  following  equation 


Vm=  f{r)pc{r)h{t  -  T)dT\t=fcs. 

J  —  OO 

Substituting  the  analog  signal  /(f)  into  (12),  we  have 

N  /*oo 

ym  =  ^2xi  /  ^i(r)pc(r)h(fcs  -  r)dr. 


(13) 


(14) 


i= l  t/-°° 

We  define  the  M  x  N  matrix  U  =  ‘b'Jf.  The  matrix  U  can  be  determined  from  the  expression  of 
ym.  It  is  clear  that  we  can  separate  out  the  expression  for  each  element  urrut  E  U  for  the  mth  row 
and  ith  column  as 

/»oo 

(15) 


uq 


/oo 

'fl’i{r)pc,rn(r)h(fcs  -  r)dr. 

-OO 


Through  the  analog-to-information  process,  rather  than  measuring  and  sampling  /(f)  directly, 
we  acquire  MciV  linear  projections  by  y  =  <fif  =  <hvI/x.  We  comment  that  the  minimal  number 
M  of  measurements  required  to  recover  the  signal  is  bounded  by 


M  >  Cp2(^}^)Klog2N 


(16) 


where  C  denotes  a  small  known  constant  (empirically  ~  2  to  2.5),  N  is  the  dimension  of  the  signal 
to  be  reconstructed,  and  /r(-,  •)  represents  the  measurement  of  the  coherence  between  <E>  and  \l>  and 
it  defined  as 

\ V)  —  max  |(0fc,  '0»)|-  (17) 

Here,  /.(2(<h,  T')  can  range  from  1  to  N .  The  smaller  /.(($.  VI/  )  is,  the  fewer  number  of  measure¬ 
ments  is  required  to  reconstruct  the  signal.  Note  that  K  represents  the  sparsity  of  the  signal,  and 
M  is  proportional  to  the  required  number  of  measurements. 

When  y  and  U  are  known,  we  have  a  common  and  practical  approach  to  find  the  sparse  solution 
to  the  following  convex  optimization  problem  [3]: 


x  =  argminxeR„  Hxlli  subject  to  y  =  Ux. 


(18) 
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where  ||  •  ||i  denotes  the  /i-norm.  The  recovered  signal  is  f  =  T'x*  with  x*  being  the  optimal 
solution  to  (18). 

We  comment  that  there  are  many  different  approaches  in  finding  the  solution  to  the  undetermined 
CS  system  described  by  y  =  Ux  =  <h\l/x.  One  common  algorithm  to  solve  the  minimization 
problem  in  (18)  is  basis  pursuit  (BP).  Since  this  algorithm  has  been  shown  to  have  a  tractable 
complexity  [3],  [15],  it  is  widely  used  in  practice. 

C.  Applications  of  CS 

CS  theory  asserts  that  one  can  recover  certain  signals  from  far  fewer  samples  or  measurements 
than  traditional  methods  at  Nyquist  rate.  In  this  section,  some  possible  applications  of  CS  is  briefly 
summarized. 

A)  Data  compression.  When  the  sparse  basis  is  unknown  or  impractical  to  implement  at  data 
compression  stage,  one  can  apply  a  randomly  designed  $  as  a  universal  encoding  strategy.  Note 
that  is  required  only  at  the  recovery  sate  of  f .  This  encoding  strategy  is  particularly  useful  for 
distributed  data  source  coding  such  as  sensor  networks  [16],  [5], 

B)  Channel  coding.  When  the  data  is  sparsity,  CS  theory  can  be  used  for  channel  coding,  e.g, 
error-correcting  code  (ECC).  Such  a  coding  system  model  can  be  described  as  [6] 

yc  =  Afc  +  e  (19) 

where  fc  is  an  input  vector,  yc  is  corrupted  measurement,  A  is  a  (coding)  matrix,  and  e  denotes 
some  unknown  error.  Based  on  the  sparsity  property,  with  proper  designing  of  A,  one  may  decode 
fc  from  measurement  yc. 

C)  Inverse  problem.  Related  discussion  is  omitted  in  this  report.  More  details  can  be  found  in 
[9]  and  [10], 

D)  Data  acquisition.  When  high  rate  or  huge  amount  of  N  DT  samples  is  difficult,  sub-Nyquist 
rate  samples  are  desired  in  designing  physical  sampling  devices  [11]. 
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E)  A/I.  It  can  be  considered  as  a  specific  application  in  data  acquisition.  As  described  in  Section 
III-B,  A/I  is  to  measure  and  compress  high  bandwidth  CT  signal  (A-sparse)  to  DT  signal  (with  a 
small  number  of  measurements). 

IV.  Through  Wall  Imaging  Radar 
A.  Preliminary  on  Through-Wall  Imaging 

Urban  radar  has  become  a  keen  research  topic  in  recent  years  because  of  the  urgent  needs  of 
urban  warfare,  counterterrorism,  calamity  rescue,  etc.  Through-wall  imaging  (TWI)  radar  systems 
can  provide  high  range  and  cross-range  solutions  by  using  ultra-wideband  (UWB)  signals  and 
antenna  arrays  with  large  apertures,  respectively. 

To  begin,  we  present  a  brief  description  of  a  conventional  TWI  radar  system.  We  first  pose  some 
assumptions: 

1.  A  target  consists  of  (at  least)  one  point  of  scatters,  and  the  reflection  coefficient  of  each  scatter 
is  a  (complex)  constant; 

2.  In  the  range  of  resolution,  human  bodies  are  replaced  by  point  targets. 

To  model  the  TWI  system,  a  stepped-frequency  scheme  is  considered.  Suppose  the  stepped 
signal  G  [fa,  /#]  has  N  narrow  band  signals  with  carrier  frequencies 

fn  —  h  +  nAf,  n  =  1,  -  ■  ■  ,  N  (20) 

where  fL  is  the  starting  frequency,  fH  is  the  highest  frequency,  and  A /  =  (  fH  —  fi)/{N  —  1)  is 
the  frequency  step. 

Figure  5  shows  the  scene  with  the  wall  for  the  point  target  V{xp,  yp)  and  the  m-th  antenna  pair 
(Tx,  /?.,.).  Suppose  that  the  wall  thickness  d  and  the  relative  permittivity  of  the  wall  are  known. 
Under  the  two  assumptions  and  the  stepped-frequency  scheme,  the  received  signal  of  frequency 
fn  at  the  receiver  position  of  the  mth  antenna  pair  can  be  represented  through  complex  envelope 


11 


Fig.  5.  A  geometry  scene  with  the  wall  for  one  target. 


(after  normalizing  by  the  transmitted  signal)  as 

p- 1 

z[m,  n]  =  ^2ap  exp(-j2nfnTpm)  (21) 

p= 0 

where  ap  is  the  reflection  coefficient  of  the  pth  target,  P  is  the  number  of  targets  in  the  scene,  m 
and  n  represent  the  spatial  index  for  m  =  0, M  —  1  and  frequency  index  for  a  =  0, N  —  1, 
respectively,  and  rpm  is  the  round-trip  propagation  delay  between  the  m-th  antenna  pair  and  the 
p- th  target.  The  round-trip  propagation  delay  can  be  calculated  through 


Tpm 


Pin.  1  “1“  -fi'ra, 2  ^  1  +  L  m,  2 


(22) 


where  Rm  j  and  Rm^  are  the  respective  distances  between  the  pth  target  and  the  refraction  points 
A  and  B,  c  denotes  the  velocity  of  light,  and  v  =  c/  s/p,-  is  the  velocity  inside  the  wall  with  sr 
denoting  the  relative  permittivity  of  the  wall. 
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We  comment  that  the  coordinates  of  the  refraction  points  A  and  B  can  be  obtained  as 


and 


Xm,l  ■I'tm 


Xm. 2  Xrm  T 


d^Xp  Xtm) 

\/£r(yp  +  d) 
d(xp  xrro) 


(23) 


(24) 


Lm, 2  xrm  i  — /  ,  ,\  • 

\J^r\Vp  +  “) 

In  the  conventional  TWI,  the  image  can  be  recovered  through  delay-and-sum  beamforming  (DS- 


BF).  Thus,  an  image  pixel  s[k,  l ]  located  at  (x/c,  yt)  is  recovered  by  the  following  DSBF  process: 

M—l  N-l 

MN 


1  M-lN-l 

2;[m,  n]  exp(j27r/nrpm) 


(25) 


m=0  n=0 


where  k  =  0, 1,  •  •  •  ,  K  —  1  and  l  =  0, 1,  •  •  •  ,  L  —  1.  In  (25),  Ff  and  L  are  the  numbers  of  pixels 
in  the  image.  Multiple- signal  classification  subspace  method  can  be  used  to  facilitate  the  image 
reconstruction  [4], 


B.  Compressed  Sensing  in  TWI  Radar 

In  TWI  radar,  if  the  number  of  targets  is  very  small  compared  to  the  number  of  pixels  in  the 
imagine,  i.e.,  P  -C  KL,  we  may  use  the  CS  technique  to  recover  such  a  sparse  target  scene. 

a)  Formulating  a  dictionary  for  TWI  Radar  Data:  First  of  all,  we  define  two  long  vector 
in  the  X  —  Y  plane  and  frequency-spatial  plane  as 

s  =  [sn,  S21,  •  •  •  ,  Sri,  Si2,  ■  ■  ■  )  SK2,  '  '  -  >  sKl\kLx  1  (26) 

and 

z  =  [pi,  321)  -  -  -  )  %li  Sl2)  '  '  -  >  SM2,  '  '  '  )  sMn]'mNx  1  (27) 

where  KL  determines  the  resolution,  s  is  a  weighted  indicator  vector,  i.e.,  the  kl  index  is  nonzero 
if  a  target  is  located  at  position  [k,l].  Here,  zmn  represents  the  data  set  z[m,n}. 

Thus,  z  can  be  expressed  as 

z  =  SFs  (28) 
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where  \P  is  a  MN  x  KL  matrix  with  \P  =  [ip1,  ip2,  •  •  •  ,  ipi,  •  •  •  ,  i^kl]-  The  wth  column  of  ^  is 


0;)«  =  exp(— j27r/aTj;U) 


(29) 


where 

\j (x,  -  a:w)2  +  2/e  +  \/(®?  -  +  2/f  ,  V7 (^w  _  x«>)2  +  d2  +  V (x&2  -  2Vb)2  +  d2) 

7i,«  = - t - 

c  c 

(30) 

where  a  =  imodiV,  b  =  \i/M\,  q  =  jmo&K,  and  e  =  [j / L\ .  Here,  |_-J  denotes  the  floor 
function,  /a  =  /z,  +  aA/,  and  A),i  and  A"/)2  can  be  found  via  (23)  and  (24).  Note  that  (xtb,  —d)  and 
(; xrb ,  0)  are  coordinates  of  the  6th  antenna  pair.  Therefore,  once  the  size  of  detect  range  and  radar 
parameters  are  specified,  the  dictionary  'h  can  be  constructed  through  (29). 

b)  Designing  of  Basis  <h:  In  CS  based  TWI  radar  system,  is  the  conversion  matrix,  while 
<b  is  the  data  acquisition  matrix.  Rather  than  sampling  z[m,n\  (or  the  MN  x  1  vector  z),  we 
measure  the  linear  projection  of  z  onto  a  set  of  basis  <h  as 


yjxl  —  MNxKLSKLxl-  (31) 

A  proper  selection  of  $  is  needed  to  minimize  the  mutual  coherence  between  <E>  and  \k.  There  are 
two  possible  selections  of  <b  are  presented. 

Principle  of  Selection  I:  Each  row  has  only  one  nonzero  element,  being  equals  to  one,  in  the 
MN  elements.  An  example  of  such  a  selection  is  shown  in  Fig.  2.  In  this  case,  both  spatial  and 
frequency  points  are  randomly  selected,  and  this  may  lead  to  relative  complex  implementation  in 
hardware  of  the  TWI  radar. 

Principle  of  Selection  II:  As  shown  in  Fig.  3,  one  first  randomly  picks  R  points  out  of  M 
antenna  pairs,  then  randomly  selects  T  out  of  N  frequency  bins  for  all  picked  R  positions  of  radar 
antennas.  This  approach  is  easier  to  implement  in  the  TWI  radar. 
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Fig.  6.  Data  acquisition  algorithm  for  the  TWI  radar  based  on  CS  (Selection  I). 
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Fig.  7.  Data  acquisition  algorithm  for  the  TWI  radar  based  on  CS  (Selection  II). 
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C.  TWl  Detection  with  CS 


Under  the  aforementioned  assumptions,  the  target  space  indicator  s  vector  can  be  recovered  from 
J  >  6'/r(<h,  \P)(log (KL/P))P  measurements.  Using  the  convex  optimization  approach,  we  can 
find  s  by  solving  an  1 1  -norm  minimization  problem 

s  =  argminj  |s|  |ix  subject  to  y  =  r&iPs.  (32) 


Since  the  dictionary  rfi,  measurement  y ,  and  the  indicator  s  are  complex,  to  facilitate  the  processing 
of  (32),  one  can  make  some  conversions  as  follows. 

Define  ©  =  rfiSP.  y  can  be  rewritten  as 


y  =  0s  =  Sft{y}  +  j%{  y} 


(33) 


where  K{y}  =  K{0}K{s}  -  3{0}3{s}  and  9f{ y}  =  3£{0}S{s}  +  3{©}3£{s}. 
Thus,  eqn.  (32)  becomes 

I  =  argmin|  |I|  |;1  subject  to  y  =  TI 


(34) 


where  y  =  [^{y}T3{y}T];[Jxl, 


1  = 


[5?{s}T3{s} 


TIT 

\2KLxl 


and 


/  »{©}  -3{©}  \ 
\  a{©}  »{©}  / 


(35) 


which  is  a  2  J  x  2 KL  matrix.  We  comment  that  as  soon  as  I  is  reconstructed,  one  can  readily 
obtain  s. 


D.  Discussion 

Based  on  our  study,  we  find  that  the  topics  on  a  CSR  system  can  be  summarized  as 

1)  Application  and  setting; 

2)  Formulation  in  the  context  of  CS:  involving  “reduced  measurement  +  reduced  DSP  complex¬ 
ity”  or  “reduced  DSP  complexity”  only; 
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3)  Design  of  the  matrix  3>; 

4)  Reconstruction  of  signal  or  target  position  (via  convex  optimization); 

5)  Impacts  from  noise  and  inaccurate  (modified)  matrix 

V.  A  General  Model  of  CS  Using  UWB 

A.  Preliminary  on  UWB 

UWB  is  one  of  the  major  breakthroughs  in  area  of  wireless  communications.  As  UWB  signal  is 
highly  sparse  in  the  time  domain,  it  is  natural  to  introduce  CS  in  the  design  of  UWB  systems. 
Typical  UWB  systems  have  the  following  features: 

1.  The  devices  have  low  power  consumption  [12]; 

2.  Signals  with  short  pulse  width  (on  the  order  of  nanoseconds). 

Based  on  the  second  feature,  one  can  find  that  the  corresponding  Nyquist  rate  could  be  on  the 
order  of  GHz.  Analog-to-digital  converters  (ADCs)  working  at  such  high  frequencies  would  be 
expensive.  Two  methods  are  proposed  to  reduce  the  sampling  frequency.  One  way  is  to  use  an 
array  of  samplers  with  different  time  offsets  to  realize  a  serial  sampling;  however,  the  hardware  cost 
will  increase  with  the  increasing  of  the  number  of  samplers.  The  other  way  is  that,  in  commercial 
applications,  UWB  receivers  may  ignore  the  details  of  UWB  signals,  i.e.,  detecting  only  the  signal 
energy  and  realizing  the  on-off  keying  (OOK)  modulation.  In  this  case,  timing  information  from 
the  peak  of  a  UWB  signal  is  lost  and  can  degrade  the  locationing  precision. 

B.  Mathematic  Model  of  UWB 

Considering  a  wireless  multi-path  channel,  the  channel  response  can  be  expressed  as 

L-l 

h(t)  =  hiS(t  -  ti )  (36) 

1=0 

where  L  is  the  total  number  of  paths,  5f)  denotes  the  impulse  function,  hi  is  the  channel  gain  of 
the  (th  path,  and  n  is  the  delay  in  path  l. 
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C.  Compression  of  UWB 

One  considers  N  +  L  —  1  transmitted  sequences  and  we  denote  the  corresponding  symbol  by 
s=  (si,  s2,  •  •  •  ,  sjv,  •  •  •  ,  sjv+l-i).  Thus,  the  sampled  output  r[nk]  at  time  n k  can  be  written  as 

L—l 

r[nk\  =  ^2  hi[nk]s[nk  -  l }  (37) 

i=o 

where  hi[nk]  is  the  /th  complex  channel  filter  tap  at  time  nk.  Assuming  that  the  channel  taps  do  not 
vary  over  N  symbol  times,  the  received  symbol  reduces  to 

L—l 

r[nk ]  =  ^2  his\nk  ~  (]  (38) 

1=0 

where  s[nk ]  =  0  for  nk  <  1.  To  analyze  the  UWB  channel,  we  transform  the  frequency-selective 
channel  into  flat  fading  channels.  As  a  result,  the  received  signal  can  be  written  as 


r  =  Hs 


(39) 


where  r  is  the  received  signal  vector  after  sampling  and  H  denotes  a  N  x  N  +  L  —  1  equivalent 
channel  matrix  which  is  given  as 


H 


hi  h0  ■■■ 
fi-2  h\  ho 


hL-i 


hi  ^ 

h*2 

h3 


(40) 


y  hL_ 1  • • •  hi  ho  J 

It  is  noted  that  (39)  is  a  standard  format  of  CS.  However,  the  matrix  H  may  not  satisfy  the  condi¬ 
tions  for  CS  signal  reconstruction  because  it  is  completely  determined  by  the  channel.  As  a  result, 
we  need  to  add  additional  compression  procedure  for  the  UWB  signal.  Such  a  compression  can  be 
done  at  either  receiver  or  transmitter  side.  Details  are  given  in  the  following  subsections. 
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c)  Receiver  Side  Compression:  Receiver  side  compression  for  UWB  signal  can  be  imple¬ 
mented  based  on  distributed  amplifier  (DA)  [13].  The  basic  operation  is  letting  the  received  UWB 
signal  pass  through  a  microchip  circuit,  along  which  the  signal  is  added  to  predetermined  weights 
and  then  sampled. 

DA,  also  called  transverse  filter,  consists  of  multiple  repeated  taps,  each  containing  a  section  of 
microchip  input  and  output  transmission  lines,  and  then  gain  cells.  Such  an  architecture  for  DA  is 
given  in  Fig.  8.  Thus,  the  output  of  a  DA  at  time  t  is 

N 

y(t)  =  cikx(t  —  hr)  (41) 

k= 1 

where  N  is  the  number  of  taps  in  a  DA,  a /3s  are  attenuation  coefficients  of  different  gain  cells, 
which  can  be  set  randomly,  x  is  the  input  signal,  and  r  is  the  fixed  time  delay  of  each  section  of 
transmission  line.  We  comment  that  the  attenuation  coefficients  and  time  delay  can  be  predeter¬ 
mined  from  hardware  design.  DA  is  suitable  for  UWB  signal  as  the  characteristic  impedance  of 
a  transmission  line  changes  very  little  over  several  GHz,  and  the  time  delays  can  easily  achieve  a 
time  scale  of  50  ps  or  less  without  changing  the  DA  structure. 

Based  on  the  proposed  DA  structure,  as  shown  in  Fig.  8,  the  received  UWB  signal  is  put  into 
the  compressor  with  M  DAs,  each  of  which  has  N  gain  cells,  and  then  sampled  by  M  ADCs.  For 
simplicity,  here  we  ignore  the  noise.  Then  the  output  of  the  ith  ADC  is  given  by 

N 

Vi  ^  ^  (fiikXk  (42) 

k= 1 

where  Xk  =  x{xTs  —  kr)  with  Ts  denoting  the  sampling  period  of  ADC  and  is  the  coefficient 
of  the  /cth  gain  cell  of  the  ith  DA. 

Let  us  define  x=  (xi,  x2,  •  ■  •  ,  xn)t  and  y=  (yi,  y2,  •  •  •  ,  Then  we  can  obtain  the  linear 

compression  equation  as 

y  =  $r  (43) 
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Fig.  8.  An  illustration  of  the  DA-based  receiver  structure. 


where  the  measurement  matrix  $  is  given  by 

/ 


$ 


011 

012 

017V 

021 

022 

027V 

0M1 

0M2  '  ' 

'  0M7V 

\ 


(44) 
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With  proper  design  of  DA,  one  can  make  (43)  suitable  for  CS  signal  processing. 

d)  Transmitter  Side  Compression:  Alternatively,  we  can  modify  the  transmitter  to  accom¬ 
modate  CS  technique.  One  approach  with  the  functionality  of  signal  mixing  is  illustrated  in  Fig. 
9. 


In  such  a  structure,  the  UWB  impulses  are  generated  according  to  the  bit  sequence  and  then  pass 
through  a  finite  impulse  response  (FIR)  filter.  The  filter  output  is  then  sent  out  through  the  wireless 
channel  and  reaches  the  receiver.  Suppose  that  the  UWB  pulse,  together  with  the  silent  period,  is 
denoted  by  an  iV-vector  s  in  the  discrete  time  domain.  Then,  the  received  signal  at  the  receiver  is 
given  by 

y0  =  f  *  h  *  s  (45) 
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Fig.  9.  An  illustration  of  the  UWB  transmitter  structure  with  signal  mixing. 


where  f  is  the  discrete  time  impulse  response  of  the  FIR  filter  at  the  transmitter  and  *  denotes  the 
convolution  operation.  Due  to  the  limitation  of  sampling  rate,  the  vector  y0  may  not  be  obtained 
directly.  After  downsampling  to  M  uniform  measurements,  the  signal  after  sampling  is  now  given 
by 

y  =  $t  *  x  (46) 

where  we  ignore  the  quantization  error,  for  simplicity  of  analysis.  The  measurement  matrix  for 
transmitter  compression  is  determined  by 

B—l 

yin)  =  X] x  inVN/M\  +  3)  (f  *  h)(-B  -  j )  (47) 

j=0 
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where  B  is  the  dispersion  of  the  composite  filter  f  *  h.  The  matrix  T»t  is  a  quasi-Toeplitz  matrix 
in  which  the  non- zero  elements  are  located  within  a  strip  and  each  column  has  B  elements.  Thus, 
with  designing  of  larger  B,  the  measurement  matrix  can  be  filled  with  more  non- zero  elements  and 
this  make  y  more  dispersed  in  the  time  domain.  More  details  of  this  approach  can  be  found  in  [14]. 


D.  Reconstruction  of  UWB 


In  UWB  system,  the  UWB  impulse  arrive  in  a  block  manner.  Hence,  in  the  unknown  vector  r,  the 
non-zero  elements  are  grouped  in  blocks,  it  is  impossible  to  have  the  non-zero  elements  distributed 
in  arbitrary  locations  of  x.  Traditional  BP  or  OMP  do  not  consider  this  prior  information.  Hence, 
by  using  block  reconstruction,  one  can  improve  the  receiver  performance. 

e)  Block  Property:  One  can  exploit  the  block  property  of  UWB  signals.  For  simplicity,  we 
ignore  the  noise  in  this  discussion.  One  can  define  the  block  sparsity  by  rewriting  the  unknown 
received  vector  r  as 


x  = 


(xi,x2,x3,--- 


(48) 


Each  subvector  x,  is  called  a  block.  Except  the  last  one,  each  xz  has  d  elements.  Here,  one 
introduces  the  block  sparsity  definition:  if  at  most  k  blocks  have  non-zero  elements,  we  say  that 
the  vector  x  is  block  /  -sparse. 

f)  li  and  l2  combined  optimization,  L-OPT:  Define  Z  =  [^],  we  denote  the  mixed  l\jl2- 
norm  as 

z 


|2  ,Z  — 


(49) 


2  =  1 

where  ||  •  ||  denotes  the  l2  norm.  The  L-OPT  minimizes  a  combination  of  l\-  and  /2-norm  of  the 
unknown  vector  x  through 

z 


min 
X  2=1 


\xz 


s.t. 


y  =  $r. 


(50) 


We  comment  that:  when  d  =  1,  one  minimizes  /2-norm  of  x;  when  d  =  N,  one  minimizes  /i-norm 
of  x.  For  a  general  value  of  d,  one  minimizes  something  between  Zi-norm  and  /2-norm. 
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g)  Condition  of  Perfect  Reconstruction:  A  matrix  A  is  said  to  satisfy  the  block  restricted 
isometric  property  (RIP)  condition  if  there  exist  parameters  5k  such  that,  for  every  vector  c  that  is 
block  A:-sparse,  such  that 

(!  -  4)||c||!  <  || Ac||l  <  (1  +  4)||c||l.  (51) 

This  condition  can  be  used  to  evaluate  the  matrix  A  in  CS-based  UWB  system  design. 

VI.  CS  Perfect  Signal  Reconstruction 

Clearly,  two  basic  conditions  are  required  for  a  CS  system.  Theoretical  analysis  of  CS  requires 
signal  sparsity.  In  addition,  RIP  or  incoherence  condition  on  sensing  matrix  is  desired  for  reliable 
signal  recovery  [17]. 

To  facilitate  the  description  of  incoherence  conditions,  we  rewrite  the  CS  signal  model  as 

y  =  <bx  =  =  ©s  (52) 

© 

where  $  is  the  measurement  matrix,  is  the  basis  matrix,  0  =  <f>,P  denotes  the  sensing  matrix, 
and  vector  s  is  A-sparse. 

A.  Conditions  on  RIP 

Let  us  define  RIP.  Suppose  that  there  exists  a  constant  5  such  that,  for  every  A -sparse  vector  \k 

(i  -  5)||v*:||!  <  ||evjr||!  <  (i  +  <j)||va-||1.  (53) 

To  continue,  we  define  the  minimum  of  all  constants  (5  satisfying  (53)  as  the  isometry  constant  5K- 
In  [19],  Davenport  and  Wakin  showed  (in  2010)  that  for 

5k  <  A=  (54) 
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if  RIP  is  satisfied,  the  orthogonal  matching  pursuit  (OMP)  algorithm  is  guaranteed  to  recover  s 
exactly  from  the  /i-norm  minimization  problem.  Later  in  2012,  Wang  and  Shim  loosened  the 


inequity  condition  in  (54)  to  [20] 


SK  < 


1 

y/K  +  1 


(55) 


B.  Conditions  on  Mutual  Coherence 

The  mutual  coherence  of  a  sensing  matrix  ©  is  defined  as 


/i(0)  =  max  |  <  0i,6j  > 


(56) 


where  6i  and  0t  are  two  properly  normalized  column  vectors  of  ©. 

Lemma  1 

If  the  following  condition  is  satisfied 


K  <  0.5 


1  + 


1 

Me). 


(57) 


OMP  algorithm  is  guaranteed  to  find  s  exactly  from  the  l  \  -norm  minimization  problem.  Note 
that  the  condition  in  Lemma  1  is  a  sufficient  and  necessary  condition.  Its  necessity  can  be  shown 
by  slightly  relaxing  the  condition  ( K  =  0.5[l  +  l//i(©)]),  which  leads  to  imperfect  recovery. 
Detailed  proofs  can  be  found  in  [21]. 

Sometimes,  perfect  recovery  may  be  hard.  Thus,  one  can  ask  for  recovery  with  high  probability. 
Such  a  condition  is  given  as  follows  [23]. 

Lemma  2 

Given  /i(@)  =  1  /y/M,  suppose  a  /.’-sparse  vector  s  follows  K  <  M/(161og(M/e))  for  some 
sufficient  small  e.  For  a  CS  system,  basis  pursuit  (BP)  algorithm  will  recover  s  with  probability 
greater  than  1  —  2e2  —  K~°  for  some  t)  >  1  s.t.  \J d  log  M /  log (M/ e)  <  c  where  c  is  an  absolute 
constant. 
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C.  Study  of  the  Coherence  Condition 


Suppose  we  have  a  CS  based  high  resolution  radar  system.  The  received  signal  can  be  formu¬ 
lated  in  the  CS  matrix  format  as 

y  =  <bx  =  <b\b  s  =  ©s  (58) 

0 

where  $  is  a  Mm  x  N  measurement  matrix,  *b  is  a  Ar  x  ML  basis  matrix,  ©  =  <b'b  denotes  a 
Mm  x  ML  sensing  matrix,  and  vector  s  is  /v -sparse. 

The  mutual  coherence  of  a  sensing  matrix  ©  is  defined  as 

/i(©)  =  max  |  <  0, .  Oj  >  |  (59) 

*7 C 

where  and  are  two  properly  normalized  column  vectors  of  ©.  One  of  our  target  is  to  find  the 
CDF  of  /i(@)  for  a  given  CS  system  (fixed  Sb). 

We  have  an  Mt  x  Mr  MIMO  CS  system  setup,  which  can  be  extended  to  an  electronic  surveil¬ 
lance  multiple-receiver  at  a  later  time.  By  discretizing  the  range-Doppler  plane  into  M  x  L  grids, 

we  are  able  to  write  the  elements  of  \b  as 

Mt 

f>uv  =  Wi(uTs )  exp (y—j2nfcTig(xv))  exp(j27rffq(uTs  -  Tiq(xv)))  (60) 

i=i 

where  1T)(-)  is  the  received  signal  waveform,  Tiq  { ■ )  denotes  the  time  delay  of  target  reflected  signal 
from  the  fih  transmitter  to  <yth  receiver  for  a  specific  grid,  and  fiq  is  the  Doppler  frequency  shift 
corresponding  to  the  grid. 

Based  on  the  CS  model,  we  derive  the  CDF  of  mutual  coherence.  We  can  rewrite  it  as  the 
maximum  of  =  \0f6j\. 

Case  I:  If  we  assume  the  elements  of  $  are  independently  drawn  from  the  random  sequence 
“{±1}”,  it  can  be  shown  that 

\gij\  =  Mm\VinI*j\.  (61) 

From  (61),  we  conclude  that  mutual  coherence  will  not  be  affected  by  <b  matrix  and  depend  on  \b 
matrix,  i.e.,  the  system  parameter. 
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Case  II:  If  we  assume  the  elements  of  $  are  independently  drawn  from  a  zero-mean  gaussian 
distribution,  the  direct  derivation  for  the  target  CDF  can  be  difficult  as  g^’ s  are  dependent  com¬ 
plex  random  variables  (RVs).  For  simplicity,  we  consider  two  elements  gi2  and  /y13  for  analytical 
analysis  purpose.  Through  a  pre-whitening  process  (detailed  steps  omitted  here),  we  have 

IO12I  =  \Tng12  +  T\2g\3 1  (62) 

and 

IO13I  =  |^2i<7i2  +  T22gi3\.  (63) 

Thus,  |^i2 1  and  |f/i3  are  found  to  be  independent  Rician  RVs. 

If  IO12I  and  |^i3 1  are  i.  n.  i.d.  RVs,  the  CDF  and  PDF  of  max{ \gi2\,  IO13I}  are  respectively  given 
by 

II  Fids)  (64) 

*=1,J=1,2 

and 

2  2 

91M  n  F-g  lk(g).  (65) 

3= 1  k=l,kd=j 

We  would  like  to  determine  a  reasonable  tight  bounding  for  the  CDF  of  max{|<7y|},  which  can 
be  used  as  a  metric  to  assess  the  recovery  performance  of  a  CS  system.  We  also  target  to  find  the 
minimum  number  of  samples  required  for  a  perfect  recovery  through  our  studies.  Such  a  parameter 
can  be  useful  in  determining  the  required  power  and  speed  of  signal  processing  unit. 

In  addition,  coherence  between  $  and  SP  may  be  investigated  to  provide  more  insights  into  the 
low  coherence  condition  for  high-resolution  CS  systems.  For  the  CS  system  discussed,  we  will 
explore  the  CS  system  with  more  realistic  conditions,  for  instance,  the  additive  noise/interference, 
in  the  future  work. 
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